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Abstract. The thick- or fat-center- vortices model has been applied to a calculation of the potentials between 
static sources of various SU (4) representations. For intermediate distances, a linear potential is obtained. 
For this region the string tensions agree qualitatively with both flux tube counting and Casimir scaling, 
even though for some representations it favors flux tube counting more. In addition, our results confirm 
the existence of two different string tensions for non-zero 4-ality representations at large distances. In this 
area, zero 4-ality representations are screened. In our computations, we have used only the first non-trivial 
vortex of 517(4). 
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' 1 Introduction 

' It has almost been proved that QCD is the true theory 
, of the strong interactions. For the high energy regime 
■ (asymptotically free region), the perturbative formulation 
shows impressive agreement between theory and experi- 
. ment. Lots of efforts have been made for the low energy 
' region where a non-perturbative formulation is required 
, and people are interested to see confinement which is one 
' of the main features of QCD. It is shown by many nu- 
merical measurements in SU{2), SU{3) and SU{4:) lat- 
tice gauge theory [1,2] that at intermediate distances a 
linear potential between the quarks for the fundamental 
and higher representations exists and the string tension is 
representation dependent and roughly proportional to the 
eigenvalue of the quadratic Casimir operator of the rep- 
resentation. This proportionality of the potential to the 
Casimir operator is called "Casimir scaling" . On the other 
hand some physicists argue for flux counting and describe 
the linear behavior of the potentials at intermediate dis- 
tances based on this idea [3,4] which claims that the string 
tension at intermediate distances is proportional to the 
number of fundamental flux tubes embedded into the rep- 
resentation. Flux tube counting is supposed to coincide 
with Casimir scaling in the large N limit. 

Besides numerical calculations, people have been try- 
ing to introduce phenomenological models to explain con- 
flncmcnt in QCD. One of these models is the center- vortex 
theory introduced in the late 1970's [5]. This model has 
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been developed by Faber et al. [6] to the thick-center- 
vortices model to study the linearity of the potential for 
higher representations and has been tested for quarks in 
SU{2). In this paper, we apply the thick- or fat-center- 
vortices model to the sources of SU{4:). The represen- 
tations 4 (fundamental), 6, 10, 15 (adjoint), 20 and 35 
of SU{4:) are studied. Each representation is shown by 
(n, m), where n and m are the number of original quarks 
and antiquarks (in the fundamental representation) antici- 
pated in producing the representation. Then since a quark 
in representation 10, for example, is constructed from two 
quarks and no antiquark: 

4(g) 4 =10 ©6. (1) 

10 is shown by (2,0). Other representations are shown as 
follows: 4: (1,0), 6: (2,0), 15: (1,1), 20: (3,0) and 35: (4,0). 
Our results from the thick-center-vortices model indicate 
that quarks in all representations of SU{A) are confined at 
intermediate distances. For large distances their behavior 
is based on the 4-ality. Representations with zero 4-ality 
are screened and other potentials will be parallel to the 
potential of either the fundamental representation or rep- 
resentation 6. I recall that for the S't/(4) gauge group, at 
large distances two different string tensions are expected. 
The general behavior of quarks in SU (4) is in agreement 
with previous works by Faber et al. [6] and Deldar [7] 
for SU{2) and SU{3), respectively. At large distances the 
potential between two sources is large enough and a pair 
of adjoint sources releases from the vacuum. Therefore 
based on the iV-ality of the representation, we expect to 
see screening or change of the slope of the potential to the 
slope of the potential in the fundamental representation. 
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This behavior has been observed for SU{2), for example, 
and reported in [6] , where at large distances the screening 
occurs for the adjoint representation and the slope of the 
potential for j = 3/2 becomes equal to that of the fimda- 
mental representation. For SU{3), the representations 3, 
6, 8, 10, 15s, 15a and 27 are studied. Results obtained for 
SU{2) and SU{3) at intermediate distances show a qual- 
itative agreement with Casimir scaling. In this paper, we 
show that for SU{4) quarks, one can get two asymptotic 
string tensions by choosing the parameters of the model 
appropriately. As studied for SU (3) [7] , the approximate 
Casimir scaling would be obtained, if one chooses a physi- 
cal profile for the vortex. For SU{A), we use one of the pre- 
tested profiles which seems physical and has worked well 
for SU{2) and SU{i). There are two different vortices for 
SU{4). For simplicity in the computations, we use only one 
of them and show that one can still get linear potentials 
for all representations at intermediate distances. Compar- 
ing the results of the SU{3) and 5J7(4) gauge groups, we 
show that even though the potentials are proportional to 
Casimir scaling, a tendency to flux tube counting is also 
observed. 

In addition, we take a closer look at the SU (3) lattice 
data for the inter-quark potentials at intermediate dis- 
tances. We discuss the; possibility of string tensions to be 
proportional to Casimir operators as well as the number 
of fundamental fluxes. 

In Sect. 2, we review very briefly the thick-center- 
vortices model, and then the results of applying this model 
to SU{4:) are given in Sect. 3. 



2 Calculation of the potentials in 5E/(4) 
by the thick-center-vortices model 

Vortex condensation theory [5] claims that the QCD vac- 
uum is flUed with closed magnetic vortices that have the 
topology of tubes or surfaces of finite thickness which 
carry magnetic flux quantized in the element of the cen- 
ter of the gauge group. In order for the vortex to have 
a finite energy per unit length, the gauge potential at 
large transverse distances must be a pure gauge. However, 
the gauge transformation which produces that potential is 
non-trivial. It is discontinuous by an element of the gauge 
center. The non-trivial nature of the gauge transformation 
forces the vortex core to have non-zero energy and makes 
the vortex topologically stable. A center vortex linked to a 
Wilson loop, in the fundamental representation of SU{N), 
has the effect of multiplying the Wilson loop by an element 
of the gauge group center, i.e. 

W(C)^exp^W{C) n = l,2,...,N-l. (2) 

The vortex theory states that the area law for the Wil- 
son loop is due to the quantum fluctuation in the num- 
ber of center vortices linking the loop. The QCD vac- 
uum does not tolerate a linear potential between adjoint 
quarks over an infinite range. Since adjoint color charges 
are screened by gluons, the force between these sources 



drops to zero and this is exactly what happens in the cen- 
ter vortex theory. However for intermediate distances, a 
linear regime is reported by lattice calculations for quarks 
of higher representations. The string tension for quarks 
of a higher representation is obtained if the vortex thick- 
ness is quite large: on the order of the typical diameters of 
low-lying hadrons. These vortices are called thick-center 
vortices [6]. The thick-center- vortices model uses two ba- 
sic assumptions to get confinement for sources of higher 
representations. The first assumption describes that the 
effect of creating a center vortex piercing the minimal area 
of a Wilson loop may be represented by the insertion of 
a unitary matrix at some point along the loop. This ma- 
trix depends on the flux distribution of the vortex and the 
generators of the group in each representation. The second 
assumption says that the probabilities / that loops in the 
minimal area are pierced by vortices and also the random 
group orientations associated with the gauge group ele- 
ments in each representation are uncorrected and should 
be averaged. Then the average Wilson loop is 

{W{C)) = n (l - E /"(I - Reer[ag(x)])| , (3) 

X [. n=l ) 

X is the location of the center of the vortex and C indicates 
the Wilson loop and Qr is defined as 

gr[cx] = ^Trexp[ia.if], (4) 

where dr is the dimension of the representation and 
{Hi, i = 1,2, ...,N — 1} are generators spanning the Car- 
tan subalgebra. / is the probability that any given unit is 
pierced by a vortex. The parameter ac{x) describes the 
vortex fiux distribution and depends on the vortex loca- 
tion; in other words, it depends on what fraction of the 
vortex core is enclosed by the Wilson loop. Therefore, it 
depends on the shape of the loop and the position of the 
center of the vortex in the plane of loop C relative to the 
perimeter. 

On the other hand, the potential may be found by 
measuring the Wilson loop and looking for the area law 
fall-off at large T: 

VF(i?,t)~exp-^(-^)^, (5) 

where R is the spatial separation of the quarks, T is the 
propagation time, and V{R) is the gauge field energy as- 
sociated with the static quark-antiquark source. For the 
thick-center-vortices model, T is assumed to be fixed and 
very huge compared to R. Therefore the loop C is just 
characterized by the width R. Thus from (3) and (5) the 
inter-quark potential induced by the vortices is 

y(i?)=^ln|l- ^/„(1-Ree.[a2r(x)])|. (6) 

X \ n— 1 ) 

V{R) depends on ac{x) (its shape and size) which is de- 
termined by the fraction of the vortex flux that is enclosed 
by the Wilson loop. 
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Vortices of type n and N — n are the same, except that 
the magnetic fluxes are in opposite directions: 

fn = fN-n and gr[aUx)]=gt[a^-''{x)]. (7) 

There arc three types of vortices in SU{4:). Because of (7), 
/i = /S) and KeGriot-ci^)] — ReSr[«c(^)] ^^'^ therefore 

{W{C)) = n {1 - 2/i(l - ReGAcxhi^)]) 

X 

-/2(1-Rea.[a2,(x)])}. (8) 

To make the computation simpler, we assume /2 = 0. This 
means that mainly the first non-trivial center element con- 
tributes and in fact it agrees with other studies of vortices 
[8] . Gr for vortex number 1 is 



Gr[oc\ = ^Trexp[i(a}iJi +a^if2 +4i?3)]- 



(9) 



The upper index of a shows the vortex profile for vortex 

number 1 which is associated with /i. If one wants to 
include vortex number 2 in the computations, one may 
choose another vortex profile with upper index 2. Vortex 
number 2 is associated with /2. The diagonal elements of 
the Cartan subalgebra in SU{A) are 



H 



1 



' ■ 2^/6 



(1,1,1,-3), 



7/2 : ^(1,1,-2,0), 
Hi : 1(1,-1,0,0). 



(10) 



Only Hi is a 4 * 4 matrix with all diagonal elements non- 
zero, and the other two come from the SU(2) and SU{^) 
gauge groups. Because of this fact, one may use only Hi 
and a\ in the calculations and thus (9) reduces to 



Gr[a\ = ^Trexp[ia}i/i]. 



(11) 



The real part of the trace of Gr[oc] in (11), with appro- 
priate normalization factor of a, contains the real part 
of the trace of (9). So again, for simplicity wc use (11). 
Calculating the normalization factor for a's is discussed 
later. 

The appropriate flux distribution, a(.T). can be chosen 
such that a well behaved potential is obtained. This means 
that with a good choice, one can see the linear term of the 
potential for all representations. In general any physical 
axially symmetric density distribution for the vortex leads 
to potentials which are acceptable in QCD. A variety of 
different fluxes is introduced in [7]. Here we choose ac{x) 
to be 

«h(x) = (^)[1 - tanh(ay(:r) + A)], (12) 

a, h and / are constants and are free parameters of the 
model. They must be chosen such that a linear potential 



at intermediate distances as well as proportionality with 
Casimir scaling are observed. In this paper a = 0.05, 6 = 4, 
/i = 0.1 and 



y{x) = 



X- R for < 
—X for \R — x\ > \x\ 



(13) 



The normalization factor for a is obtained based on the 
assumptions of the model. Equation (2) shows that each 
time that a vortex links to a fundamental Wilson loop, 
the loop is multiplied by the factor exp^r^ . Since we are 
using only one vortex and N = 4, we have: 



W{C) exp* W{C). 



(14) 



Thus, based on the thick-center-vortices model, for higher 
representations, every time that the minimal surface is 
pierced by a center vortex, a center element exp(^) should 
be inserted somewhere along the loop: 



W{C) = Tr[UU...U] 



rr[i7C/...exp[y]...;7]. 



(15) 



Using this fact and the point that we have used only one 
vortex and only Hi in our computations, the normaliza- 
tion factor is obtained from exp^ = exp'"^'^ which hap- 
pens if the vortex core is entirely contained within the 
Wilson loop. Therefore a should satisfy the following con- 
ditions. 

(1) Vortices which pierce the plane far outside the loop do 
not affect the loop. In other words, for fixed i?, as a; — s- oo, 

(2) If the vortex core is entirely contained within the loop, 
then a = V^ir. 

(3) As i? ^ then a ^ 0. 

Equation (12) with its normalization factor satisfies 
the above conditions. Changing the parameters a and b 
by any factor F and setting a ^ and b ^ bF only 
changes the scale of the potential and the physics remains 
the same. However not any arbitrary value for a and b 
leads to a linear potential. 

Figure 1 shows this flux distribution versus x, the lo- 
cation of the center of the vortex. For this plot R = 50. 
Thus the size of the loop is such that the vortex may 
overlap completely the Wilson loop. Another example of 
a physical flux distribution is the one which is zero every- 
where except on the boundary of the vortex [7] . A density 
proportional to a delta function which is zero everywhere 
except at the two points where the vortex first enters and 
exits the Wilson loop is among the non-axially symmet- 
ric distributions. Using this flux, one loses Casimir scaling 
proportionality. 



3 Results and discussion 

Figure 2 shows the potentials for various representations 
versus R in the range R E [1, 100]. For all representations, 
there is a region where the potentials are linear (about 
R e [5,15]). In general, at large distances, screening for 
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Fig. 1. Flux distribution of (12). For this plot i? = 50 and 
the Wilson loop may entirely overlap the vortex core 



15 (adjoint) and 35 should occur. 4-ality which is {n, m) 
mod 4, is zero for representations 15: (1, 1) and 35: (4,0). 
Screening is expected for zero 4-ality representations. In 
fact, at large distances, the potential between two sources 
is large enough to create a pair of quarks in the adjoint 
representation (15) from the vacuum. Then the original 
static sources combine with the 15s and for zero 4-ality 
representations a 15 component is created and the poten- 
tial will be screened. For non-zero 4-ality representations, 
two possible scenarios may occur: for example, for rep- 
resentation 20, a fundamental component will be created 
and therefore the slope of the potential changes to that 
of the fundamental representation; on the other hand, for 
6 and 10 dimensional diquark representations, a quark in 
representation 6 is created and the slope of the potential 
in representation 10 changes to the slope of the potential 
in the representation 6. We have 

15 o 6 = 64 e 10 e 10 e 6, (16) 
15 10 = 7006401096. (17) 

This is why there are two different string tensions, one for 
the 4 and another for the 6 dimensional representation. 
For the SU{N) gauge group with N > 3, there is more 
than one asymptotic string tension. In S'f7(4), potentials of 
the 6 and 10 dimensional representations do not parallel 
the fundamental representation potential and the slope 
of the potential of representation 10 changes to that of 
representation 6 [10]. Using only vortex number 1 in our 



Fig. 2. Static sources potential for the range of i? £ [1, 100]. 

At intermediate distances, potentials arc linear. At large dis- 
tances, representations 15 (adjoint) and 35 are screened; rep- 
resentation 20 gets the same slope as the fundamental repre- 
sentation and representation 10 is parallel to representation 
6 



Table 1. This tabic shows Casimir numbers ratios, number 
of flux tubes and string tensions ratios for the SU{4) gauge 
group. A qualitative agreement of string tensions with both 
Casimir scaling and flux tube counting is observed 



Repn. 4 (fund.) 


6 15a 


10 20 35 


(n,m) (1,0) 


(2,0) (1,1) 


(2,0) (3,0) (4,0) 


Cr/Cf 1 


1.33 2.13 


2.4 4.2 6.4 


fund, fluxes 1 


2 2 


2 3 4 


kr/kf 1 


1.51 1.56 


1.76 2.31 2.66 



calculations, /2 = 0, we have gotten two different string 
tensions for the 5?7(4) sources. 

Figure 3 plots the ratios of the potential of each repre- 
sentation to that of the fundamental one. Although these 
ratios start up roughly at the ratios of the correspond- 
ing Casimirs, which are 1.3, 2.4, 2.13, 4.2 and 6.4 for the 
representations 6, 10, 15, 20 and 35, respectively, but in 
the most linear part of the potential (Fig. 4) , the ratios of 
some representations are closer to the flux tube counting. 
Table 1 shows the representations we have studied with 
the number of original quarks and antiquarks, {n,m), an- 
ticipated in each representation. In the third row, the ratio 
of Casimir number of each representation to that of the 
fundamental one is indicated. The number of fundamen- 
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Table 2. The same as Table 1 but for the SU (3) gauge group 



Repn. 



3 (fund.) 8 6 15a 10 27 15s 



(n,m) (1,0) 

Cr/Cf 1 

fund, fluxes 1 

kr/kf 1 



(1,1) (2,0) (2,1) (3,0) (2,2) (4,0) 
2.25 2.5 4 4.5 6. 7 
2 2 3 3 4 4 

2.02 2.21 3.1 3.4 3.8 5.6 



to 



o 



_ 1 1 1 1 1 1 1 

'-,35/! 


II 


1 1 1 _ 


20/f 










10/f 




6/f 






II 


II 


II 


II 



10 

R 



15 



20 




Fig. 4. The most linear part of the potentials. The slope of 
each potential is obtained for the range R G [7, 12]. They are 
.282, .427, .439, 498, .652 and .752 for representations funda- 
mental, 6, adjoint, 10, 20 and 35, respectively. The ratios of 
string tensions (slopes of the linear part of potentials) are re- 
ported in Table 1 



Fig. 3. Potential ratio of each representation to the fundamen- 
tal representation. These ratios start up roughly at the ratios 
of the corresponding Casimirs, which are 1.3, 2.4, 2.13, 4.2 and 
6.4 for representations 6, 10, 15, 20 and 35, respectively, but 
in the most linear part of the potential (Fig. 4), the ratios are 
also in agreement with the flux tube counting especially for the 
higher representations. The slope of the potentials will become 
constant after the confinement regime is finished. It happens 
at about ii = 30 with the given parameters (o, b, /) we used in 
the model 



tal fluxes existing in each representation is shown in the 
fourth row. The last row indicates the ratio of the string 
tension of each representation to the fundamental one, cal- 
culated from this work. These string tensions are obtained 
from the range of R G [5, 12] which is the most linear part 
of the potentials in our calculations. The agreement with 
both Casimir scaling and flux counting is qualitative, but 
it favors flux tube counting more. Comparing Tables 1 and 
2, one can sec that the tendency to the flux tube count- 
ing increases by increasing the number of gauge groups, 
especially for the representations 35 of SU (4) and 15s of 
SU{3). 

Looking closer at Fig. 3, it seems that there arc four 
regions for the potentials. For the first area, basically for 



R < 5, the potentials are proportional to Casimir scaling; 
especially for very small R. The potential ratios start out 
at the ratios of Casimirs which are larger than the number 
of fundamental fiuxes, but they get close to the number of 
fltixes embedded in each representation at about R = 10 
for the representations 35, 20 and 10 and at about i? = 20 
for representations 15 (adjoint) and 6. Based on the as- 
sumptions of the thick-ccnter-vortices model, vortices arc 
uncorrelated and also the orientation of the "vortex inser- 
tion" in color space is chosen at random. Thus, it might 
be possible that these assumptions lead effectively to the 
interaction of flux tubes. Then one may explain the; bc;- 
havior of the potentials at different distances by the in- 
teraction between fundamental strings as follows: When 
the distance between sources is very small, the fundamen- 
tal fluxes overlap and the string tension is larger than 
the number of fluxes times the fundamental string ten- 
sion. This is because of the repulsion between the fluxes 
[4]. This may explain the behavior of the potential for 
distances less than 5 in Fig. 3. On the other hand as the 
distance increases, the fluxes tend to attract each other 
and therefore a negative energy is added to the binding 
energy of fluxes and this makes the string tension smaller 
(for 10 < i? < 20) [4,9]. In general if there is no interaction 
between the fluxes, the string tension of the representa- 
tion must be Kaf where K is the number of fundamental 
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4 Conclusion 



10 



X string tension ratios from lattice 
o Casimir ratios 



No. of flux tubes 



Representation 



Fig. 5. Ratios of string tensions of 517(3) quarks of different 
representations to the string tension of quarks in the funda- 
mental representation are plotted. Considering the lattice data 
errors, good agreement with both Casimir scaling and flux tube 
counting is observed 



Using thick-center-vortices model, we have shown that for 

all representations at intermediate distances there is a re- 
gion where quarks are confined and the potentials are lin- 
ear. For this linear regime, our results show evidence of 
proportionality of the string tensions with both Casimir 
scaling and flux counting even though the agreement is 
slightly better for the flux tube counting. At large dis- 
tances, zero 4-ality representations (15, 35) are screened, 
the potential of the representation 20 has become parallel 
to that of the fundamental representation and the poten- 
tials of the diquark representations (6 and 10) get the 
same slope. In our computations, we have used only the 
first non-trivial vortex of SU{4). Although we have not 
used vortex number 2 in our calculations, our results at 
intermediate distances agree with lattice results and the 
potentials are linear and qualitatively in agreement with 
Casimir scaling at intermediate distances. 

We are doing new computations with /2 not necessarily 
zero and study the effect of vortex number 2 and also the 
different core sizes of the vortices in the linear region of 
the potential. Especially, we are interested in studying the 
effect of different parameters of the model on removing 
the concavity of the potentials which has been observed 
for SU (2) and SU (3) sources, as well. 



fluxes and a j is the fundamental string tension. Screening 
or change of the slope of the potentials starts from about 
i? = 30. 

Casimir scaling obtained from SU{3) lattice calcula- 
tions by Bali [3] has been verified for distances less than 
Ifm and the ratios are reported to be bigger than the 
number of fluxes. It seems that the thermal distance be- 
tween strings in the lattice results is not large enough and 
an overlap between the strings leads to a repulsion force 
and therefore increasing of the string tension. This coin- 
cides with the iZ < 5 of Fig. 2. It should be noticed that 
the scales of R and V(R) are arbitrary (adjustable) in 
the thick-center-vortices model. On the other hand, the 
results reported by Deldar [7] are obtained for larger dis- 
tances (up to 2.5 fm). Figure 5 is plotted based on the 
data of [7]. Cross signs show the ratios of the string ten- 
sions of SU (3) sources of various representations to that of 
the fundamental representation. Circles indicate Casimir 
ratios, and diamonds show the number of fundamental 
strings in each representation. As claimed by the lattice 
people, the ratios of the string tensions are proportional to 
the Casimir operators. On the other hand, the plot shows 
that there is a rough agreement with the number of funda- 
mental tubes as well. One may conclude that the errors of 
the lattice data are still too large to discriminate between 
the two hypotheses, Casimir scaling and flux counting. 



5 Acknowledgement 

We would like to thank Dr. S. Olejnik and Dr. M. Faber 
for their valuable comments and patience in answering our 
questions. We are grateful to Dr. C. Bernard for his help, 
especially his suggestion of studying SU{4:) potentials. We 
also thank the research council of the University of Tehran 
for support of this work. 



References 

1. C. Bernard, Phys. Lett. B 108, 431 (1982); Nucl. Phys. 
B 219, 341 (1983); J. Ambj0rn, P. Olcsen, C. Peterson, 
Nucl. Phys. B 240, 189 (1984); C. Michael, Nucl. Phys. B 
259, 58 (1985); G. Poulis, H. Trottier, Phys. Lett. B 400, 
358 (1997); S. Ohta, M. Wingate, Nucl. Phys. 83 (Proc. 
SuppL), 381 (2000) 

2. G. Bah, Phys. Rev. D 62, 114503 (2000); S. Deldar, Phys. 
Rev. D 62, 034509 (2000) 

3. C. Michael, Adiabatic surfaces from the lattice excited 
gluonic potentials, hep-ph/9809211; G.S. Bah, Phys. 
Rept. 343, 1 (2001) 

4. A. Armoni, M. Shifman, Nucl. Phys. B 671, 67 (2003) 

5. G 't Hooft, Nucl. Phys. B 153, 141 (1979); J.M. Cornwall, 
Nucl. Phys. B 153, 392 (1979); G. Mack, V.B. Petkova, 
Ann. Phys. (NY) 123, 442 (1979); 125, 117 (1980); Z. 
Phys. C 12, 177 (1982); H.B. Nielson, P. Olesen, Nucl. 
Phys. B 160, 380 (1979); J. Ambj0rn, P. Olesen, Nucl. 
Phys. B 170, 265 (1980); B 170, 60 (1980); L.G. Yaffe, 
Phys. Rev. D 21, 1574 (1980) 

6. M. Faber, J. Greensite, S. Olejnik, Phys. Rev. D 57, 2603 
(1998) 



S. Deldar, S. Rafibakhsh: SU{4) string tensions 



7 



7. S. Deldar, JHEP 0101, 013 (2001) 

8. L. Del Debbio, D. Diakono, Phys. Lett. B 544, 202 (2002) 

9. L. Del Debbio, M. Faber, J. Greensite, Nucl. Phys. B 414, 

594 (1994) 

10. S. Ohta, M. Wingatc, Nucl. Phys. 73 (Proc. SuppL), 
435 (1999); B. Lucini, M. Tcpcr, Phys. Lett. B 501, 128 
(2001); L. Del Debbio, H.H. Panagopoulos, P. Rassi, E. 
Vicari, JHEP 0201, 9 (2002) 



